The entanglement entropy of one-dimensional gases 
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We introduce a systematic framework to calculate the bipartite entanglement entropy of a spatial 
subsystem in a one-dimensional quantum gas which can be mapped into a noninteracting fermion 
system. To show the wide range of applicability of the proposed formalism, we use it for the 
calculation of the entanglement in the eigenstates of periodic systems, in a gas confined by boundaries 
or external potentials, in junctions of quantum wires and in a time-dependent parabolic potential. 
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Entanglement is a fundamental phenomenon of quan- 
tum mechanics. Much theoretical work has focused on 
the entanglement properties of quantum many-body sys- 
tems, showing their importance to characterize the many- 
body dynamics P^j. In particular, lots of studies have 
been devoted to quantify the highly nontrivial connec- 
tions between different parts of an extended quantum 
system, by computing von Neumann or Renyi entangle- 
ment entropies of the reduced density matrix of a 
subsystem A. The most remarkable result is the univer- 
sal behavior at ID quantum critical points, determined 
by the central charge of the underlying conformal field 
theory (CFT) [2H5]. For a partition of an infinite ID 
system into a finite piece A of length £ and the remain- 
der, the entanglement entropy for £ much larger than the 
short-distance cutoff a is 

Si = -TY[pA\npA]^^hi-+0{l), (1) 
6 a 

where c is the central charge. This result has been con- 
firmed in many spin-chains and in ID itinerant systems 
on the lattice [Tl. These studies have allowed a deeper 
understanding of ID simulation algorithms based on the 
so-called matrix product states However, the same 
result must be valid also in systems in continuous space 
(when the UV cutoff is properly imposed). Apart from 
the interest to describe trapped ID gases experimentally 
realized with cold atoms, the entanglement of continuous 
models is also instrumental to develop ID tensor network 
algorithms for g the one proposed in [7J. Despite 

of this fundamental interest, almost no effort (with the 
exception of Refs. [8j [9] and the orbital partitioning in 
quantum Hall states [10]) has been devoted to the spatial 
entanglement of gas models (that is distinguished from 
the particle partitioning |llj). 

In this Letter we present a systematic framework to 
tackle free gases in any external conditions for an arbi- 
trarily large number of particles. The most general re- 
sult of this investigation is that, when dealing with a fi- 
nite number of particles N, the ID entanglement entropy 
grows like In TV, with a prefactor that again is given by 
the central charge. In this formulation N acts as an ex- 
plicit UV cutoff, representing a concrete alternative to 



the lattice. While the IniV behavior could have been 
predicted also on the basis of scaling arguments, its va- 
lidity in many different physical situations is not obvi- 
ous a priori. Furthermore, being the method very gen- 
eral, it allows us to derive exact predictions that, in some 
cases, were not known from the lattice and/or CFT argu- 
ments. We apply our framework only to gases of spinless 
fermions, but its validity is more general. It is indeed 
straightforward to include spin degrees of freedom. Fur- 
thermore, the ID Bose gas in the limit of strong short- 
ranged repulsive interaction (corresponding to a gas of 
impenetrable bosons, describing also the dilute limit of 
the finite-strength model can be mapped to spinless 
fermions, so that their entanglement entropies coincide. 

In the following we first present the general framework 
based on a mapping from a continuous Fredholm determi- 
nant into a standard one of dimension N. We then show 
the power of the method applying it to several physical 
situations, including wires junctions and non-equilibrium 
problems that are not yet solved on the lattice. 

The method. Let us consider a system of N non- 
interacting spinless fermions with discrete one-particle 
energy spectrum. As well known, the many body wave 
functions ^'(a;i, xn) can be built from the one-particle 
eigenstates via ^'(xi, xjv) = det[(/)fc(x„)]/A/]vT, where 
the normalized wave functions (j)k{x) represent the oc- 
cupied single-particle energy levels. The ground state is 
obtained by filling the lowest N energy levels. Thus, the 
ground-state two-point correlator reads 

N 

C{x, y) ^ {cHx)c{y)} = ^ cl>l{x)Mv) , (2) 
fe=i 

where c{x) is the fermionic annihilation operator and the 
one-particle eigenfunctions (t>k{x) are ordered according 
to their energies. 

We want to compute the bipartite entanglement en- 
tropy of a space interval A, extending from xi to X2, 
in this fermion gas. For this purpose, we introduce the 
Fredholm determinant 

Da{X) =dct[AI-q , (3) 

where the continuous matrix C and the identity I are re- 
stricted to the part at hand, i.e., from xi to X2. Then, 
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calculations based on the Wick theorem [T3], and the 
integral representation of Ref. '14] for the Renyi entan- 
glement entropy, allow us to write 



In Trp^ 
1 - a 



d\ din £>a(A) 



where the integration contour encircles the segment [0, 1], 
and 



«(A) = 



1 



In [A" + (1 - A)"] 



(5) 



For a —7- 1 it reduces to the von Neumann definition ([ij. 
The integral representation Q has been already derived 
and used in the context of discrete chain models [13] , thus 
involving the determinant of a standard matrix with the 
lattice sites as indeces. 

The Fredholm determinant is turned into a standard 
one by introducing the N x N overlap matrix A (also 
considered in Ref. 81) with elements 



dz(f>l{z)(f)miz), n,m = l,...,N, (6) 



such that TrC'' = TrA*-' for any k, thus 



fc=i 



fcA'^ 



E 

k=l 



ln(A-a™) , 

rn— 1 

(7) 

where the eigenvalues of A [we use the stan- 

dard relation for Fredholm determinants In det [I — zM] = 
— X]fc^i(-^'^Tr M''')/fc and we drop off a term cx In A giv- 
ing vanishing contribution in Eq. Q]. Inserting it into 
the integral Q, we obtain 



N 



N 



Sa{xi,X2) = / ^ E \ ^}^a E (8) 

m— 1 m— 1 

as a consequence of the residue theorem. 

The matrix A is easily obtained for any non-interacting 
model from the one-particle wave functions, as the defini- 
tion ([g]) shows. Calculating the entanglement entropies 
is then reduced to an iV x iV eigenvalue problem that 
can be easily solved numerically and in some instances 
even analytically, as we are going to show. Details of the 
calculations will be reported elsewhere. 

Ground state of a periodic system. In a system of 
length L with periodic boundary conditions (BC), the 
one-particle wave-functions are (j)k{x) = e^'"^^^^^ /^/L 



with wave numbers A: S Z. In the ground state of the 
fermion gas, filling the N fc-modes with lowest energies, 
the matrix A associated with a segment of length I is 



sin7r(r7, — m)X 
Tr{n — m) ' 



X = £/L, 



(9) 



n, m = 1, N. By inserting the N eigenvalues of A into 
Eq. ^ , we obtain the entanglement entropy in a system 






N=20 
□ N=19 
N=50 
. N=100 
N=200 




1 , 



12 3 4 
In 2Nsin(jiX) 



0.1 



0.2 0.3 
X 



0.4 



FIG. 1: (a) Entanglement entropy Si for periodic (uppermost 
set of data) and Dirichlet (bottom) BC compared with the 
asymptotic results for different values of A'^ {X = £/L). (b) 
The function F2{X) for the finite-size scaling behavior of the 
particle-hole excited state vs the CFT prediction |20j . 



of N particles. Furthermore, since In Da — Indet G, with 
G = AI — AanA^xA^ Toeplitz matrix, we can use the 
Fisher-Hartwig conjecture [T3] to rigorously infer that 

Sa = 1(1 + a-^)\n(2N simrX) + ba + O(N-i), (10) 
6 

where ba and the leading 0{N~i) corrections can also 
be computed analytically [TB]. Fig.[T]shows a comparison 
with exact finite- iV calculations, for a = 1. Eq. (10) 



agrees with the CFT prediction for finite systems [4^, 
obtained from Eq. ([iJ replacing £ with the chord length 
L/tt sin(7rX). Thus, Eq. (10 1 represents the first explicit 



analytic confirmation of this CFT prediction. 

Systems with hard-wall potential. We now consider a 
gas of spinless fermions confined in the interval [0, L] by 
a hard- wall potential (i.e., the gas density vanishes for 
X ^ [0,i]). We consider a bipartition starting from the 
boundary, i.e., A = [0,i], but the general case can also 
be treated. Some algebra similar to before leads to 



sin[7r(n — m)X] sin[7r(n -I- m)X] 



7r(n — m) 



7r(n -|- to) 



(11) 

Using a recent generalization of the Fisher-Hartwig con- 
jecture to Toeplitz-f Henkel matrices [121 HZj: we obtain 

Sa = ^(l + «-i)ln(47Vsin7rX) + ^ + 0(7V-i). (12) 



A comparison of the finite-A^ results with Eq. ( |12[ ) is 
shown in Fig. [l] The leading oscillating 0{N~i) cor- 
rections can also be computed. They correspond to the 
0{£-i) corrections found within CFT {TE\\T9\. 

Excited states can be easily treated in this formalism 
by summing in Eq. (|2| over the occupied one-particle 
states representing the excited states. As an example, 
let us consider the particle-hole excitation in a periodic 
system obtained by moving one particle from the high- 
est occupied level to the first available one. The cor- 
responding N X N matrix A has the first n — 1 raws 
and columns identical to the ground state, and the last 
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different. Note that it ceases to be a Toeplitz matrix. 
In Fig. [T] we compare the scaling function F2{X) = 
exp[S'2 — l/41n(2A^sin7rX) — ^2] for several values of N 
with the corresponding CFT prediction [20] F2{X) — 
[7+cos(27rX)]2/64. The CFT curve is clearly approached 
in the large- limit. 

Star graphs or wire- junctions. A novel application of 
the outlined method consists in determining the entan- 
glement at a junction of quantum wires (called also star 
graph). Networks of quantum wires with junctions at- 
tracted recently a lot of attention |^ , mainly because of 
their possible applications in nanocircuits. A star graph 
consists of M wires joining in a single point (called ver- 
tex). We consider wires whose bulk is described by non- 
interacting spinless fermions. The only interaction is lo- 
calized at the vertex and it is encoded in a M x M scatter- 
ing matrix S [HJ , as pictorially showed in Fig. [2j We 
consider wires of finite length L with hard-wall bound- 
ary conditions at their ends. Scale invariant junctions 
can be constructed and classified [HJ [52] : they are either 
isolated points or families with free parameters. 

Let us first consider the case of two wires, i.e., a gas of 
free fermions with a localized impurity at the vertex. The 
allowed scale invariant conditions at the vertex are [22] : 
the two trivial ones (Dirichlet and Neumann) which dis- 
connect the two wires and give no entanglement, and the 
one-parameter family described by the scattering matrix 



He) 



1 



-1 + e' 
2e 



2e 



1 



(13) 



(A phase can be added, but it does not enter any mea- 
surable quantity and will be omitted.) Note that e = 1 
corresponds to full transmission, i.e., no impurity. 

In order to derive the ground-state entanglement en- 
tropy of one wire, we compute the matrix A in Eq. ([6| 
by diagonalizing the one-particle problem. We obtain 

" ,(1/2) forriT^m, (14) 



A. 



where 



1 + 
1 



1 



for odd 71, 



1 



for even n. 



is defined in Eq. (11). Then, we determine 



the wire entanglement entropy Sa{s; N) as a function of 
the transmission coefficient 

(15) 



1 + e^ 

We observe a logarithmic growth 
Scis;N) =C„(s) In TV 



Oil): 



(16) 



with a prefactor Ca{s) that depends on s and not only 
on the central charge, see Fig. [2] Assuming universality 
in terms of s, we can exploit the results of Ref. 23J for 
the lattice Ising and XX models with a defect, to infer 

1 + g-2Qw(s) 

Ca{s) = -^Tvi r / dxln' 



7r2(l 



(1 + e-2t^(s))c 



(17) 




FIG. 2: Left: M wires interacting through the scattering ma- 
trix S. Right: Entanglement entropy 5*1 for a junction of two 
wires and for different values of the transmission coefficient 
s. The full lines correspond to the asymptotic behavior 1 16 1 
with Ci(s) given by the limit a — >■ 1 of Eq. (17 1. 



where a;(s) = acosh [(coshx)/s]. The perfect agreement, 
see Fig. [2] with the large- TV behavior of our numerical 
evaluations nicely confirms the universality conjecture. 

The case of a general number of wires M can be anal- 
ogously treated. After diagonahzing the model [22], we 
construct the correlation matrix reduced to one wire a, 
and then the 2N/M x 2N/M matrix A (for simplicity we 
assume iV to be a multiple of M) corresponding to the 
same wire. After some algebra we obtain 



A„™-T„v/1-T2 B„„(l/2) 
A.„„ Tl for even n, A„„ = 



for n ^ m, (18) 
(1 - T^) for odd n. 



where = (1-|-Saa) /2. Notice that A has the same form 
of that for the two- wire problem, cf. Eq. (14 1, which is 
recovered setting Ta = e/Vl + e^- Thus, using Eq. (16 1, 



we also derive the asymptotic behavior of the entangle- 
ment entropy of any of the M wires given by Eq. ( 16 ) 



with s ^ 2Ta^l-T2 = ^l-§2^. Note that s2is 
the total probability of a signal to be transmitted from 
the wire a to all other wires. Being Ca a monotonous 
function of s, the maximum entanglement entropy is al- 
ways obtained for s = 1 (§aa = 0), i.e., in the absence of 
reflection. 

Non-equilibrium evolution. All results presented so far 
concern the eigenstates of free-fermionic Hamiltonians. 
However, the method allows us to also obtain numerically 
exact and/or analytic computations out of equilibrium. 
Indeed Eq. ^ is valid if we replace the one-particle eigen- 
functions (pkix) with appropriate solutions of the corre- 
sponding time-dependent one-particle Schrodinger equa- 
tion. As an example that can be treated analytically, we 
consider the off-equilibrium quantum dynamics of a gas 
in a time-dependent harmonic potential. 

Using the method outlined above and the knowl- 
edge of the eigenstates of the harmonic oscillator, one 
can straightforwardly calculate the ground-state (equi- 
librium) entanglement entropy of any spatial subsystem 
[a;i,X2]. In particular, the half-space entanglement en- 
tropy, i.e., of [— cx),0], behaves as S'q(— oo,0) = -^{1 + 
a~^)hiN -\- ea-\-0{N~a)^ which can be analytically ob- 
tained, including the constant Bq, by developing results 
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of Ref. [Uj . In the presence of a time- dependent poten- 
tial V{x) = ^K{t)x'^, and starting from the ground state 
for a given kq — K,{t — 0), we can exploit the solution of 
the one-particle problem [551 US] to infer that the time- 
dependent entanglement entropy behaves as 



Sa{xi,X2;t) = Saixi/s{t),X2/s{t);0) 



(19) 



where s{t) is an analytical function of the time-dependent 
potential with s(0) = 1. This shows the remarkable prop- 
erty that the evolution in a harmonic potential simply 
corresponds to a global rescaling of the system size. In 
the case of a quantum quench with an instantaneous re- 
moval of the harmonic potential of frequency ^/kq, we 
have s{t) = vT+^^o^- Then for t — )• oo, s(t) di- 
verges and the entanglement of any finite interval van- 
ishes in the long time limit, while the entanglement of 
any semi-infinite piece [— oo,a;] tends asymptotically to 
^„(-oo,0;0). 

Other non-equilibrium situations such as local quan- 
tum quenches injunctions of quantum wires (e.g., instan- 
taneously turning on/ofF the point contact at the vertex) 
can also be tackled within this framework. They may 
provide important insights in view of the recent proposals 
of using the full counting statistics after a quench as an 
experimental probe and a measure of entanglement [37] . 

Discussions. In this Letter we have introduced a gen- 
eral framework to calculate the entanglement entropy of 
gases that have a representation in terms of free fermions. 
The asymptotic entanglement entropy grows like InA'^, 
with a prefactor that is given by the central charge and 
that (in known cases) agrees with universal predictions. 
The UV regularization given by the number of particles 
A'' is alternative to the lattice spacing and, for the de- 
scription of gases, it has an immediate physical interpre- 
tation. Moreover, the formulas of the universal features 
of the entanglement are easier to handle in the continuum 
at finite N than on the lattice. 

We conclude by mentioning some possible extensions 
of this work beyond the straightforward (but still inter- 
esting) applications to other ID free gases. Motivated by 
recent experiments in cold atoms, it is interesting to con- 
sider ID gases in an external random potential that, in 
the case of a perfect gas, can be treated with our method. 
All the results we have presented are in ID, but the ap- 
plication of the method to higher dimensions is straight- 
forward. 
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